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Observation of Top Quark Production in pp Collisions 

Abstract 

We establish the existence of the top quark using a 67 pb _1 data sample 

' of pp collisions at y/s = 1.8 TeV collected with the Collider Detector at Fer- 

0^ ■ milab (CDF). Employing techniques similar to those we previously published, 

H ■ we observe a signal consistent with ti decay to WWbb, but inconsistent with 
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the background prediction by 4.8<r. Additional evidence for the top quark is 
provided by a peak in the reconstructed mass distribution. We measure the 



^ ' top quark mass to be 176 ± 8(stat.) ± 10(sys.) GeV/c 2 , and the tt production 

(N ■ 

O ' cross section to be 6.8 94 pb. 
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Recently CDF presented the first direct evidence for the top quark[|T|, the weak 
isodoublet partner of the b quark required in the Standard Model. We searched for 
tt pair production with the subsequent decay tt — »■ WbWb. The observed topology in 
such events is determined by the decay mode of the two W bosons. Dilepton events 
(e/i, ee, /i/i) are produced primarily when both W bosons decay into ev or \w. Events 
in the lepton+jets channel (e, /i+jets) occur when one W boson decays into leptons 
and the other decays into quarks. To suppress background in the lepton+jets mode, 
we identify b quarks by reconstructing secondary vertices from b decay (SVX tag) 
and by finding additional leptons from b semileptonic decay (SLT tag). In Ref. 
we found a 2.8a excess of signal over the expectation from background. The inter- 
pretation of the excess as top quark production was supported by a peak in the mass 
distribution for fully reconstructed events. Additional evidence was found in the jet 
energy distributions in lepton+jet event s0. An upper limit on the tt production 
cross section has been published by the DO collaboration 0. 

We report here on a data sample containing 19 pb _1 used in Ref. |1| and 48 pb _1 
from the current Fermilab Collider run, which began early in 1994 and is expected to 
continue until the end of 1995. 

The CDF detector consists of a magnetic spectrometer surrounded by calorime- 
ters and muon chambers[|J. A new low-noise, radiation-hard, four-layer silicon vertex 
detector, located immediately outside the beampipe, provides precise track recon- 
struction in the plane transverse to the beam and is used to identify secondary ver- 
tices from b and c quark decays 0. The momenta of charged particles are measured in 
the central tracking chamber (CTC), which is in a 1.4-T superconducting solenoidal 
magnet. Outside the CTC, electromagnetic and hadronic calorimeters cover the pseu- 
dorapidity region \r)\ < 4.20 and are used to identify jets and electron candidates. 
The calorimeters are also used to measure the missing transverse energy, E't, which 
can indicate the presence of undetected energetic neutrinos. Outside the calorimeters, 
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drift chambers in the region \r]\ < 1.0 provide muon identification. A three-level trig- 
ger selects the inclusive electron and muon events used in this analysis. To improve 
the tt detection efficiency, triggers based on Er are added to the lepton triggers used 
in Ref. 0. 

The data samples for both the dilepton and lepton+jets analyses are subsets of a 
sample of high-PT inclusive lepton events that contain an isolated electron with Ex > 
20 GeV or an isolated muon with P T > 20 GeV/c in the central region (|ry| < 1.0). 
Events which contain a second lepton candidate are removed as possible Z bosons if an 
ee or fifi invariant mass is between 75 and 105 GeV/c 2 . For the lepton+jets analysis, 
an inclusive W boson sample is made by requiring 

t T > 20 GeV. Table | classifies 
the W events by the number of jets with observed Et > 15 GeV and \r]\ < 2.0. The 
dilepton sample consists of inclusive lepton events that also have a second lepton with 
Pt > 20 GeV/c, satisfying looser lepton identification requirements. The two leptons 
must have opposite electric charge. 

The primary method for finding top quarks in the lepton+jets channel is to 
search for secondary vertices from b quark decay (SVX tagging). The vertex- finding 
efficiency is significantly larger now than previously due to an improved vertex-finding 
algorithm and the performance of the new vertex detector. The previous vertex- 
finding algorithm searched for a secondary vertex with 2 or more tracks. The new 
algorithm first searches for vertices with 3 or more tracks with looser track require- 
ments, and if that fails, searches for 2-track vertices using more stringent track and 
vertex quality criteria. The efficiency for tagging a b quark is measured in inclusive 
electron and muon samples which are enriched in b decays. The ratio of the mea- 
sured efficiency to the prediction of a detailed Monte Carlo is 0.96 ± 0.07, with good 
agreement (±2%) between the electron and muon samples. The efficiency for tagging 
at least one b quark in a tt event with > 3 jets is determined from Monte Carlo to 
be (42 ±5)% in the current run, compared to the (22 ±6)% reported in the previous 
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publication]^]. In this paper we apply the new vertex finding algorithm to the data 
from the previous and the current runs. 

In Ref. [0] , we presented two methods for estimating the background to the top 
quark signal. In method 1, the observed tag rate in inclusive jet samples is used to 
calculate the background from mistags and QCD-produced heavy quark pairs (bb and 
cc) recoiling against a W boson. This is an overestimate of the background because 
there are sources of heavy quarks in an inclusive jet sample that are not present 
in W+jet events. In method 2, the mistag rate is again measured with inclusive 
jets, while the fraction of W+jet events that are Wbb and Wcc is estimated from 
a Monte Carlo sample, using measured tagging efficiencies. In the present analysis, 
we use method 2 as the best estimate of the SVX-tag background. The improved 
performance of the new vertex detector, our ability to simulate its behavior accurately, 
and the agreement between the prediction and data in the W + 1-jet and W + 2- 
jet samples make this the natural choice. The calculated background, including the 
small contributions from non-W background, Wc production, and vector boson pair 
production, is given in Table [I]. 

The numbers of SVX tags in the 1-jet and 2-jet samples are consistent with the 
expected background plus a small ti contribution (Table |l| and Figure |l|). However 
for the W + > 3-jet signal region, 27 tags are observed compared to a predicted 
background of 6.7 ± 2.1 tags|§. The probability of the background fluctuating to 
> 27 is calculated to be 2 x 10 -5 (see Table ||) using the procedure outlined in 
reference 1 ||. The 27 tagged jets are in 21 events; the 6 events with 2 tagged jets 
can be compared with 4 expected for the top+background hypothesis and < 1 for 
background alone. Figure [l| also shows the decay lifetime distribution for the SVX 
tags in W + > 3-jet events. It is consistent with the distribution predicted for b decay 
from the ti Monte Carlo simulation. From the number of SVX tagged events, the 
estimated background, the calculated ti acceptance, and the integrated luminosity of 
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the data sample, we calculate the tt production cross section to be 6.81^4 pb, where 
the uncertainty includes both statistical and systematic effects. This differs from the 
cross section given in Ref. by 6.9 ± 5.9 pb. 

The second technique for tagging b quarks (SLT tagging) is to search for an 
additional lepton from semileptonic b decay. Electrons and muons are found by 
matching CTC tracks with electromagnetic energy clusters or tracks in the muon 
chambers. To maintain acceptance for leptons coming directly from b decay and from 
the daughter c quark, the Pt threshold is kept low (2 GeV/c). The only significant 
change to the selection algorithm compared to Ref. |l] is that the fiducial region for 
SLT muons has been increased from \r)\ < 0.6 to \r)\ < 1.0, resulting in an increase of 
the SLT total acceptance and background by a factor of 1.2. 

The major backgrounds in the SLT analysis are hadrons that are misidentified 
as leptons, and electrons from unidentified photon conversions. These rates and 
the smaller Wbb and Wcc backgrounds are determined directly from inclusive jet 
data. The remaining backgrounds are much smaller and are calculated using the 
techniques discussed in Ref. The efficiency of the algorithm is measured with 
photon conversion and J/i/j — > data. The probability of finding an additional e or 
fi in a tt event with > 3 jets is (20 ±2)%. Table ||] shows the background and number 
of observed tags for the signal region (W+ > 3 jets). There are 23 tags in 22 events, 
with 15.4 ± 2.0 tags expected from background. Six events contain both an SVX 
and SLT tag, compared to the expected 4 for top+background and 1 for background 
alone. 

The dilepton analysis is very similar to that previously reported |lj, with slight 
modifications to the lepton identification requirements to make them the same as 
those used in the single lepton analysis. The dilepton data sample, described above, 
is reduced by additional requirements on Ex and the number of jets. In order to 
suppress background from Drell-Yan lepton pairs, which have little or no true JZl T , the 
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FJ T is corrected to account for jet energy mismeasurement [|I|] . The magnitude of the 
corrected % T is required to be at least 25 GeV and, if FJ T is less than 50 GeV, the 
azimuthal angle between the FJ T vector and the nearest lepton or jet must be greater 
than 20°. Finally, all events are required to have at least two jets with observed 
E T > 10 GeV and \rj\ < 2.0. 

The major backgrounds are Drell-Yan lepton pairs, Z — > tt, hadrons misiden- 
tified as leptons, WW, and bb production. We calculate the first three from data and 
the last two with Monte Carlo simulation Jl|. As shown in Table ||] the total back- 
ground expected is 1.3 ± 0.3 events. We observe a total of 7 events, 5 e/i and 2 \x\x. 
The relative numbers are consistent with our dilepton acceptance, 60% of which is in 
the e/x channel. Although we have estimated the expected background from radiative 
Z decay to be small (0.04 event), one of the /i/i events contains an energetic photon 
with a /i/i7 invariant mass of 86 GeV/c 2 . To be conservative, we have removed that 
event from the final sample, which thus contains 6 events. Three of these events con- 
tain a total of 5 6-tags, compared with an expected 0.5 if the events are background. 
We would expect 3.6 tags if the events are from tt decay. When the requirement that 
the leptons have opposite charge is relaxed, we find one same-sign dilepton event (efi) 
that passes all the other event selection criteria. The expected number of same sign 
events is 0.5, of which 0.3 is due to background and 0.2 to ti decay. 

In summary, we find 37 6-tagged W+ > 3-jet events that contain 27 SVX tags 
compared to 6.7 ± 2.1 expected from background and 23 SLT tags with an estimated 
background of 15.4 ± 2.0. There are 6 dilepton events compared to 1.3 ± 0.3 events 
expected from background. We have taken the product (P) of the three probabilities 
in Table ^| and calculated the likelihood that a fluctuation of the background alone 
would yield a value of P no larger than that which we observe. The result is 1 x 10~ 6 , 



which is equivalent to a 4.8a deviation in a Gaussian distributionjKj. Based on the 
excess number of SVX tagged events, we expect an excess of 7.8 SLT tags and 3.5 
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dilepton events from ti production, in good agreement with the observed numbers. 

We have performed a number of checks of this analysis. A good control sample 
for 6-tagging is Z+jet events, where no top contribution is expected. We observe 15, 
3, and 2 tags (SVX and SLT) in the Z + 1-jet, 2-jet, and > 3-jet samples respectively, 
compared with the background predictions of 17.5, 4.2, and 1.5. The excess over 
background that was seen in Ref. is no longer present. In addition, there is no 
discrepancy between the measured and predicted W + 4-jet background, in contrast 
to a small deficit described in reference 1 (Tl| . 

Single lepton events with 4 or more jets can be kinematically reconstructed 
to the ti —>■ WbWb hypothesis, yielding for each event an estimate of the top quark 
mass|[]]. The lepton, neutrino (J^ T ), and the four highest-Ex jets are assumed to be the 
ti daughters |T2|. There are multiple solutions, due to both the quadratic ambiguity in 
determining the longitudinal momentum of the neutrino and the assignment of jets to 
the parent Ws and 6's. For each event, the solution with the lowest fit x 2 is chosen. 
Starting with the 203 events with > 3 jets, we require each event to have a fourth jet 
with Et > 8 GeV and < 2.4. This yields a sample of 99 events, of which 88 pass 
a loose x 2 requirement on the fit. The mass distribution for these events is shown in 
Figure |2|. The distribution is consistent with the predicted mix of approximately 30% 
ti signal and 70% W+jets background. The Monte Carlo background shape agrees 
well with that measured in a limited-statistics sample of Z+4-jet events as well as in a 
QCD sample selected to approximate non-W background. After requiring an SVX or 
SLT 6-tag, 19 of the events remain, of which 6.9^;g are expected to be background. 
For these events, only solutions in which the tagged jet is assigned to one of the b 
quarks are considered. Figure ^| shows the mass distribution for the tagged events. 
The mass distribution in the current run is very similar to that from the previous 
run. Furthermore, we have employed several mass fitting techniques which give nearly 
identical results. 
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To find the most likely top mass, we fit the mass distribution to a sum of the 
expected distributions from the W+jets background and a top quark of mass M top [|TJ. 
The -ln(likelihood) distribution from the fit is shown in the Figure^ inset. The best fit 
mass is 176 GeV/c 2 with a ±8 GeV/c 2 statistical uncertainty. We make a conservative 
extrapolation of the systematic uncertainty from our previous publication, giving 
M top = 176 ± 8 ± 10 GeV/c 2 . Further studies of systematic uncertainties are in 
progress. 

The shape of the mass peak in Figure |] provides additional evidence for top 
quark production, since the number of observed 6-tags is independent of the observed 
mass distribution. After including systematic effects in the predicted background 
shape, we find a 2 x 10~ 2 probability that the observed mass distribution is consistent 
with the background (Kolmogorov-Smirnov test). This is a conservative measure 
because it does not explicitly take into account the observed narrow mass peak. 

In conclusion, additional data confirm the top quark evidence presented in 
Ref. [|TJ. There is now a large excess in the signal that is inconsistent with the 
background prediction by 4.8<r, and a mass distribution with a 2 x 10 -2 probability 
of being consistent with the background shape. When combined, the signal size and 
mass distribution have a 3.7 x 10~ 7 probability of satisfying the background hypoth- 
esis (5. Oct). In addition, a substantial fraction of the jets in the dilepton events are 
6-tagged. This establishes the existence of the top quark. The preliminary mass and 
cross section measurements yield M top = 176 ± 8 ± 10 GeV/c 2 and a t i = 6. 8^2*4 pb. 
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Fermilab staff. We thank the staffs of our institutions for their many contributions 
to the construction of the detector. This work is supported by the U.S. Department 
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Ministry of Education, Science and Culture of Japan, the National Science Council 
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observed 


observed 


background 


Njet 


events 


SVX tags 


tags expected 


1 


6578 


40 


50 ± 12 


2 


1026 


34 


21.2 ±6.5 


3 


164 


17 


5.2 ±1.7 


> 4 


39 


10 


1.5 ±0.4 



Table 1: Number of lepton±jet events in the 67 pb _1 data sample along with the 
numbers of SVX tags observed and the estimated background. Based on the excess 
number of tags in events with > 3 jets, we expect an additional 0.5 and 5 tags from 
ti decay in the 1 and 2 jet bins respectively. 



Channel: SVX SLT Dilepton 



observed 


27 tags 


23 


tags 


6 events 


expected background 


6.7 ±2.1 


15.4 


±2.0 


1.3 ±0.3 


background probability 


2 x 10" 5 


6 x 


io- 2 


3 x 10~ 3 



Table 2: The numbers of tags or events observed in the three channels along with the 
expected background and the probability that the background would fluctuate to the 
observed number or more. 
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Figure 1: Number of events before SVX tagging (circles), number of tags observed 
(triangles), and expected number of background tags (hatched) versus jet multiplicity. 
Based on the excess number of tags in events with > 3 jets, we expect an additional 
0.5 and 5 tags from tt decay in the 1 and 2 jet bins respectively. The inset shows the 
secondary vertex proper time distribution for the 27 tagged jets in the W + > 3-jet 
data (triangles) compared to the expectation for b quark jets from tt decay. 
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Figure 2: Reconstructed mass distribution for the W + > 4-jet sample prior to 
6-tagging (solid). Also shown is the background distribution (shaded), with the nor- 
malization constrained to the calculated value. 
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Figure 3: Reconstructed mass distribution for the 6-tagged W + > 4-jet events (solid). 
Also shown are the background shape (dotted) and the sum of background plus tt 
Monte Carlo for M top = 175 GeV/c 2 (dashed), with the background constrained to the 
calculated value, 6.9tf;g events. The inset shows the likelihood fit used to determine 
the top mass. 
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